Flux induced semiconducting behavior of a quantum network 
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We show that a diamond shaped periodic network, recently proposed as a model of a spin filter 
(A. Aharony et al. [l[) is capable of behaving as a p-type or an n-type semiconductor depending 
on a suitable choice of the on-site potentials of the atoms occupying the vertices of the lattice, and 
the strength of the magnetic flux threading each plaquette of the network. A detailed study of the 
density of states of an infinite network is made together with the conductance of finite sized system 
to establish the idea. 

PACS numbers: 63.50.+X, 63.22.+m 
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Low dimensional model quantum systems have been the 
objects of intense research, both in theory and in exper- 
iments, mainly due to the fact that these simple looking 
systems are prospective candidates for nano devices in 
electronic as well as spintronic engineering p], 0, d, BL 
H B 0, B S M, HO, 51 0- Apart from this feature, 
several striking spectral properties are exhibited by such 
systems owing to the quantum interference which is spe- 
cially observed in quantum networks containing closed 
loops. Examples are, the Aharonov-Bohm (AB) effect [j| 
in the magnetoconductance of quantum dots [7( , electron 
transport in quantum dot arrays [IIJ3, Fano effect in a 
quantum ring-quantum dot system [81, spin filter effects 
in mesoscopic rings [HI [ll[ and dots @ , to name a few. 

Recently, Aharony et al. [H,0| have proposed a model of 
a nano spintronic device using a linear chain of diamond- 
like blocks of atomic sites. Each plaquette of the array is 
threaded by identical magnetic flux. They have analyzed 
how the Rashba spin orbit interaction and the AB flux 
combine to select a propagating ballistic mode. A similar 
chain was earlier investigated by Bercioux et al. [H, [l3| 
in the context of spin polarized transport of electrons. 
However, there are certain special spectral features of- 
fered by the diamond chain, particularly the role of the 
AB flux, which we believe, remain unexplored. This is 
precisely the area we wish to highlight in the present 
communication. We show that an infinite diamond chain 
of identical atoms behaves as an insulator at T = K in 
the presence of a non-zero AB flux. As we arrange atoms 
of two different kinds (represented by two different values 
of the on-site potential) periodically on a diamond chain, 
a highly degenerate localized level is created near one of 
the two sub-bands of extended states. The proximity of 
this localized level to either of the sub-bands can be con- 
trolled by tuning the AB flux, and can be made to stay 
arbitrarily close to either of the sub-bands. The entire 
system is then capable of behaving as an n-type or a p- 
type semiconductor as explained later. The conductance 
spectrum of a finite array of the diamond plaquettes is 
also studied to judge the applicability of such a network 
geometry in device engineering. 

We adopt a tight binding formalism, and incorporate 
only the nearest neighbor hopping. We begin by refer- 



ring to Fig. [lja). A magnetic (AB) flux (j> threads each 
plaquette. The Hamiltonian of the network is given by, 
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where, Cj (c|) are the annihilation (creation) operator at 
the ith site of the network, ej is the on-site potential 
at the i-th site which we shall choose as ea or es as 
shown, and t is the constant nearest neighbor hopping 




FIG. 1: (a) Schematic view of a section of an infinite diamond 
chain, (b) The RG scheme. The arrow gives the direction of 
the forward hopping texp(i8ij). 

integral. The phase dij is constant, and is given by, #y = 
±27r0/40o, where <po = hc/e is the flux quantum. The 
positive or negative sign of the phase factor depends on 
whether the electron hops from site A(B) to site B(A) 
along the arrow (called the 'forward' hopping) or against 
it ('backward' hopping). 

To obtain the average density of states (AVDOS) of the 
infinite system we make use of the system of equations 
satisfied by the Green's function Gy, viz, 
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On the right hand side of Eq. [2] the index k runs over 
the nearest neighbors of the i-th site. The density 
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of states is obtained by a decimation renormalization 
group (RG) method The RG scheme is depicted 

in Fig. [TJ The lower ^-vertices in Fig. [JJa) are deci- 
mated first, and the original diamond array gets trans- 
formed into an array of triangular plaqucttcs with the 
surviving A atoms sitting at the top and the pair of 
original £>-atoms getting renormalizcd into a pair of B\ 
atoms, each characterized by the effective on-site poten- 
tial eflj = e B + 2t 2 /(E— ca). The Bi-Bi effective hopping 
integral is Tf = t 2 exp(i7r</>/0o)/ (E — ca) and its complex 
conjugate t\, — t f . Xf = t exp(i7T0/20o) is the B±-A 
hopping integral as shown in Fig. [1Tb). R s conjugate is 
named Xb- The subscripts '/' and 'b' refer to the forward 
or backward hopping due to the broken time reversal sym- 
metry due to the magnetic flux. In Fig.[JJb) we have not 
shown the flux explicitly as the flux has been automati- 
cally included in Xf, Xb, Tf and Tf,. The array of triangles 
is now renormalizcd by decimating the alternate A sites 
and the appropriate -B-sites preserving the triangular ge- 
ometry of the array. The recursion relations satisfied by 
the parameters are given by, 



e A = e A + p 2 x b + PiXf 
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for the site potentials and, 
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for the hopping integrals between Bi-A and B\-B\ pairs 
on the renormalized triangular array. Obviously, x' b and 
r' b are given by, xf and tJ* respectively, and, pi = 



[x b {E - e 2 ) + x f Tf}/5, p 2 = p*, h f = Tf + x\j(E - e A ), 



h b = h* f , n 
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r\ with e 2 — 
e A ) and S = (E - e 2 ) 2 - t/t 6 . With 
a small imaginary part added to the energy E, the local 
Green's functions at the sites A and B of the original 
diamond array are then obtained as, 
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where, the superscript * above refer to the fixed point 
values of the respective parameters. The AVDOS is given 
by 



p(E) = ^pa(E) + ^p B (E) 



(6) 



where Pi{E) = — (l/ir)Im Gu, i being A or B. 

Let us now present two separate cases which will throw 
light on the central problem addressed in this paper, viz, 
the semiconducting behavior of such a network. 

Case I: €a = e_B 

The AVDOS in this case is illustrated for e A = eg = 
in Figs. [DJa), (b) and (c) for = 0, 4> — 0o/5 and for 
4> = 2</>o/5 respectively. In the zero flux case the spec- 
trum is a continuum with a very high value of the density 



of states at the center, i.e. at E — 0. The system ex- 
hibits a metallic character. As the magnetic field through 
the plaquette assumes non-zero value, a gap opens at the 
centre with the peak in the AVDOS fixed at E = 0. This 
high value of the AVDOS corresponds to strongly local- 
ized states. It is to be appreciated that these states are 
localized strictly on the A-type vertices of the array and 
the care has to be taken while renormalizing the diamond 
array so that both the A-type sites do not get decimated. 
This precise RG scheme is explained earlier. If the state 
corresponding to an energy E is localized, then the renor- 
malized hopping integral will iterate to zero. This is ob- 
served for the central peak in the AVDOS by iterating 
Eq. 0] The gap around the central peak widens as the 
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FIG. 2: p(E)-E for an infinite diamond network array when 
the AB flux is (a) cf> = 0, (b) 4> = 0o/5, and (c) <^ = 20 o /5. We 
have chosen ea = es = and t = 1. The AVDOS has been 
plotted within the range zero to two. 



flux is increased gradually, and the sub-bands at the two 
extremities shrink to two sharp lines of zero width as 
(\> = 4>o/2. The picture reverses as 4> increases from </>o/2 
to 4>q. The exactly symmetrical location of the local- 
ized level with respect to the sub-bands at the flanks 
makes the system behave as a semiconductor when the 
electronic filling of the system lies between 1/3 and 2/3. 
More precisely, one can get a p-type semiconductor when 
the filling factor n e is 1/3 and, an n-type semiconductor 
when 1/3 < n e < 2/3. 

Case II: / «b 

The AVDOS is exhibited in Fig.[3]with e B = and e A = 2 
for different values of the flux. By comparing Fig. [2] and 
Fig. [3] we see that the highly degenerate localized level 
is pinned at E — 6a- This has been verified using var- 
ious values of keeping eb fixed. So, the states are 
localized at the vertices A of the diamond chain. The 
other interesting feature is that as the gaps open up for 
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non-zero flux, the localized level is placed asymmetrically 
with respect to the continuous sub-bands of extended 
eigenstates. With a proper tuning of the flux, one such 
sub-band can be brought arbitrarily close to the sharp 



The strongly degenerate level is pinned at E — 2. If the 
magnetic field is small, then the energy gap between the 
localized level at E = 2 and the bottom of the conduction 
band (the right sub-band) is small enough for electrons to 
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FIG. 5: Schematic view of a finite diamond chain attached to 
two metallic electrodes. 

bridge. The system now behaves as an n-type semicon- 
ductor. By reverting to the case of tA = — 2 we can sim- 
ulate a p-type semiconductor-like behavior observed by 
the same diamond chain. In this case the localized level 
is pinned at E = —2 and, we need to fix the Fermi level 
at the top of the valence band, so that electrons can jump 
into unoccupied levels above creating holes in the valence 
band. It should be noted however, that with increasing 



FIG. 3: p(E)-E for an infinite diamond network array when 
the AB flux is (a) <j> = 0, (b) <j> = cj> /5, and (c) 0=20 o /5. We 
have chosen €a = 2 and e_g = with t = 1. The AVDOS has 
been plotted within the range zero to two. 



localized level. For example, with €a — 2, the localized 
level resides closer to the upper sub-band (the conduction 
band) than the lower one (the valence band). A rever- 
sal of the sign of the on-site potential €a reverses the 
picture. In either case, the localized level can be placed 
arbitrarily close to any one of the sub-bands by appropri- 



2 




FIG. 4: Variation of the electron concentration n as a function 
of the flux <f). The red, green and blue curves correspond to 
the temperatures ksT — 0.5 , 1 and 1.5 respectively, fcs being 
the Boltzmann constant. 

ately tuning the value of the magnetic flux. This control 
over the proximity of the degenerate localized level and 
a sub-band can be utilized to simulate an extrinsic semi- 
conductor like behavior. For example, let us refer to the 
case in which tA = 2 and es — 0, and fix the Fermi 
level at E — 2 (i.e. all states from the bottom of the left 
(valence) band up to E — 2 are filled up at T = K). 
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FIG. 6: g-E (red color) and p-E (blue color) curves for a 20- 
ring chain, (a) cj> = 0, (b) cf> = 0.2, and (c) (j> = 0.4. Other 
parameters are, ea = es = 0, t = 3, eo = 0, to = 4 and 
ts = td = 2.5. 

flux the gap-width increases and hence the probability 
of electrons crossing over to any one of the sub-bands 
decreases. 

In Fig. 2] we present the variation of the concentra- 
tion of electrons (holes) at finite temperatures in the con- 
duction (valence) bands for €a = 2 and —2 respectively. 
With flux increasing from zero the carrier concentration 
(electron or hole) diminishes. 

The applicability of the above physics depends on 
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whether such features are exhibited by systems with a 
finite size as well. In order to check this we have calcu- 
lated the density of states and conductance of a twenty 
(20) plaquette chain. To calculate the conductivity, the 
system is connected between two metallic electrodes, viz, 




— 8—4 o 4 8 

E 




— 8 — 4 O 4 8 




— 8 — 4 O 4 

E 



FIG. 7: g-E (red color) and p-E (blue color) curves for a 20- 
ring chain, (a) <j> = 0, (b) <j> = 0-2, and (c) <j> = 0.4. Other 
parameters are, = 2, es = 0, t = 3, eo = 0, to = 4 and 
TS — TD = 2.5. 

the source and the drain (Fig. [5]), described by the stan- 
dard tight binding Hamiltonian and parametrized by a 
constant site potential eo and nearest neighbor hopping 
integral t Q . The hopping integral between the source 
and the system is rg, while it is T£> between the system 
and the drain. Throughout the calculation we choose 
the units where c = e = h = 1. For low bias and 



low temperature one calculates the conductivity using 
the single channel Landauer formula [l5| g = (2e 2 /h)T, 
where, the transmission coefficient T is given by [151 ] 
T = Tr[T s G r r D G a }. T s and r D correspond to the 
imaginary parts of the self energy due to the coupling 
of the diamond chain with the electrodes, and G repre- 
sents the usual Green's function. In Figs. [5] and [7J we 
have shown simultaneously the AVDOS and the conduc- 
tance spectra. Fig. [6] illustrates the symmetric case with 
ca = es = 0. With zero flux, the system exhibits oscil- 
lating conductance profile with resonance peaks. As the 
flux becomes non-zero, and increases in value, the con- 
ductance windows occupy the regions corresponding to 
the two sub-bands at the flanks, and shrink in width as 
the flux increases towards the half flux quantum. Fig- 
ure [7] depicts a similar qualitative feature, but now with 
an asymmetric conductance profile as the potentials ca 
and €b are different. 

Before we end, it may be mentioned that whether an 
array of such plaquettes will behave as an insulator or 
will have a metallic character in the absence of an exter- 
nal magnetic field, depends on the geometry of the unit 
cells. For example, with an array of identical triangular 
plaquettes a gap already exists in the middle of the spec- 
trum even when the external magnetic field is zero. Thus 
such a triangular array will have an insulating character 
at T = K and at zero flux. However, by changing the 
site potentials one can generate strictly localized levels 
inside the gap and can again have the n- or p-type semi- 
conducting behavior as discussed earlier in respect of the 
diamond plaquette-array. 

Finally, in view of the potential application of the net- 
works as a device we would like to point out that the 
qualitative features presented here should remain valid 
even at finite temperature (~ 300 K) since, the broad- 
ening of the energy levels of the diamond array due to 
its coupling with the electrodes will be much larger than 
that of the thermal broadening [lj| . 
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